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1. Homogeneity of the cytosol approximation: comparison of shell 
models 

One major assumption in our modeling approach is the homogeneity approximation of the cytosol. That 
is, we assume that the proteins of interest are evenly distributed within this compartment, and that 
diffusion is infinitely fast. This is a valid approximation if the diffusion rate is much faster than the 
observed transport rates. If this were not the case, models that account for a more complex structure of 
the cytosol should give a better fit to the data. In this analysis we have tested if models in which we have 
subdivided the cytosol into several smaller compartments give significantly better fit than our simple 
two-compartment model. 

We have approached this issue by constructing models where the cytosol is divided into a series of 
concentric shells. The nucleus is the innermost shell; it is encompassed by the innermost shell of the 
cytosol, which is in turn encompassed by the next layer of the cytosol, and so on. In this analysis we have 
compared dividing the cytosol from 1 (the standard two-compartment model) to 3 such shells. By 
definition each shell of the cytosol has been constructed to have equal volume (and therefore slightly 
different thickness). 

The models have a diffusion constant, which in the models corresponds to the transport rate between 
neighboring shells. This constant is higher between outer shells than inner shells, in proportion to the 
relative differences in surface area contacting the shells. The constant in the model is a free parameter, 
constrained to be larger than 1 and smaller than 10, and it's a multiplication factor of the first diffusion 
parameter. The optimization algorithm has chosen values between 1 and 8, depending on the data set. 
For the data set we show here was ~6. 
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We thus have the following models: 

• 1-shell: 10 parameters (2 kinetic, 8 initial conditions)* 

• 2-shell: 15 parameters (3 kinetic, 12 initial conditions) 

• 3-shell: 20 parameters (4 kinetic, 16 initial conditions) 

*ln the main paper we have actually used a reduced version with 6 parameters, of which only four 
are initial conditions. 

For each shell added, we add 1 diffusion parameter, and 4 initial conditions parameters. This 
purposefully gives the model a lot of freedom, for example, we allow for it to completely change the 
state values between FRAP-trains in an unrealistic and discontinuous way. This is because if we can then 
reject these versions of the model, we can also reject versions where more care has been taken into 
choosing realistic initial conditions of the various shells. 

We analyzed the three models on ten different data sets in the following way. We classified our ~90 data 
sets into 5 different categories, based on the observed kinetics. These were: slow dynamics, medium 
dynamics, fast dynamics, linear response, and no response. From each category we chose two 
representative data sets, one with relatively low x 2 cost, and one with relatively high x 2 cost.We fitted the 
more complex models to the same data sets and observed the improvement in fit. We then tested if this 
improvement was significant using a standard likelihood ratio test. The results can be seen in the table 
below. 



Class 


X 2 Cost 1-shell 


Decrease 2 -shell 


Decrease 3 -shell 


Slow - high 


147.62 


-6.62 


-2.55 


Slow - low 


47.67 


-6.25 


-0.85 


Medium - low 


71.31 


-1.26 


-1.10 


Medium - high 


113.16 


-2.49 


-1.92 


Fast - low 


107.59 


-2.48 


-6.83 


Fast - high 


176.89 


-18.10 


-0.91 


Nothing - high 


317.26 


-8.13 


-0.53 


Nothing - low 


151.07 


-1.79 


-2.82 


Linear - low 


52.81 


-2.09 


-1.55 


Linear - high 


97.71 


-7.09 


-5.04 


Cutoff 




11.07 


11.07 



Table SMI. Shown are x 2 costs for the ten different data sets. First column shows the actual cost for the 
simplest model. The second and third columns show the successive decrease in x 2 cost gained by adding 
a new shell. This value should be compared to the cutoff-value for a likelihood ratio test, which follows a 
X 2 distribution. Shown is the cutoff for 5% significance level. Significant results are highlighted in orange. 
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From our statistical analysis, we could only see a significant improvement in one of the cases (see Table 
1). Note that with a 5% significance level we expect to see one false positive one time out of twenty, 
which is precisely the frequency we observe above. However, we choose to study the observed data set 
in more detail to see if we can understand how the more complex model achieves this better fit. The 
results can be seen in Figure SMI. The more complex model does not fit to the data in a qualitatively 
different way, rather it accomplishes the improve cost by matching the initial conditions more flexibly. 

In conclusion, one in ten data sets showed an improvement in goodness of fit when going from a 1-shell 
model of the cytosol to a two-shell model. However, for this data set, the improvement did not seem to 
come from a qualitatively different fit to the data. No data set showed any significant improvement 
when going from a 2-shell to a 3-shell model. Taken all together, it does not seem that modeling the 
cytosol in a more complex way improves the overall goodness of fit. 
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Figure SMI. Comparison between a 1-shell (top) and 2-shell (bottom) model of the cytosol for the data 
set that showed an improved goodness of fit. Red error bars show the observed cytosolic data, blue 
error bars show the observed data in the nucleus. Dark red and dark blue continuous lines show model 
fits to this data for the 1-shell and 2-shell model top and bottom respectively. Continuous green and 
magenta lines in bottom figure show the time course for the inner and outer shell respectively in the 2- 
shell model. 
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2. Theoretical analysis on the impact of noise in the accuracy of the 
determined parameters 



First we simulated data similar to the one we get with our experimental protocol. To do this, we used the 
same ODEs 2-compartment model than for the fittings, with the time resolution of the experiments and 
"realistic" parameters and initial conditions (kl=0.4, kEV=4.41, n 0 = 500, c 0 =500). Then we introduced 
some noise according to the following equation 

simData 

noisyData = simData + r * 

snr 

Where snr represents the signal to noise ratio (that was set to 1, 5, 10 , 50 or 100) and r is a random 
number generated with the randn function of MATLAB (that returns random numbers, normally 
distributed with mean=0 and sd=l). 

For each signal to noise ratio, we generated 10 datasets (Figure S2). When we studied the result of 
varying the number of FRAPs per train, we generated 10 independent datasets for each train length from 
1 to 10 (Figure la). Then, the simulated data was treated the same way than experimental data. 

3. FRAP quantification and correction 

We corrected the obtained data for the non-cellular fluorescence background, cells autofluorescence, 
imaging photobleaching (FigSM2) and reversion of the intentional photobleaching (Fig SM3). 

Imaging photobleaching 

We tested two approaches to reduce the photobleaching: one is to perform partial FRAPs, and another is 
to reduce the number of images acquired after the bleaching. We analyzed both of these scenarios with 
simulated data and concluded that they don't have a significant impact on the estimation of the 
parameters (not shown). Based on all this analysis, we decided to perform trains of four partial FRAPs 
(bleaching approximately 50% of the signal), and to image the recovery for only 7 sec. Also, since we 
don't need a small confocal section, and in order to use the lowest possible laser power during imaging, 
we set the confocal microscope pinhole to wide open (500 u,m). Even under these conditions, the 
imaging photobleaching was still detectable and we had to quantify it to correct the data (Figure SM2). 
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Figure SM2. Quantification of the photobleaching due to the imaging of the cells. Left, curves of average 
fluorescence over time for different laser intensities. Other than the laser intensity, the conditions were 
set as in the FRAP experiments (5 pictures per second, pinhole wide open). Right, characteristic time of 
decay (tau) for the curves on the left as a function of the laser intensity. 

Reversibility of photobleaching 

Essential for a FRAP experiment is the assumption that any increase in fluorescence signal be due to 
movement of molecules to the bleached region and not from other processes, such as fluorophore 
maturation and reversibility of the photobleaching process. Maturation of YFP is too slow to have a 
significant effect during our FRAPs (half time is ~ 40 min (Gordon et al., 2007)). However, it was 
important to determine and correct for the fraction, if any, of the photobleached YFP that could revert to 
the fluorescent state (Sinnecker et al., 2005; Henderson et al., 2007). 

The reversibility of photobleaching was measured by completely photobleaching whole cells (so that 
fluorescence could not recover by moving from other compartments) using laser intensities between 10 
and 100% (Figures SM3 a and b). For these experiments, cells were pre-treated with cycloheximide for at 
least for 3 hours before the experiments, to allow all fluorescent molecules to mature (Gordon et al., 
2007). We found that laser intensity did not change the characteristic time of the recovery, but it 
strongly changed the fraction of fluorescence that recovered (Figures SM3 c and d). 

Remarkably, for the YFP used here, 20% of the bleached molecules recovers its fluorescence after 1 min 
when photobleached with a 100% laser intensity as we used for FRAP experiments (figure SM3 c). In our 
experiments of trains of 4 FRAPs, the time lapse between the first photobleaching pulse and the last 
imaged time-point is 27 sec, and the fluorescence recovered in that period corresponds to approximately 
10% of the bleached molecules. 

Similar experiments using cells fixed with 4% PFA showed no recovery (not shown). The discrepancy of 
these results can be explained if the fixating process also affects the recovery of the photobleaching and 
thus suggest that this control cannot be performed on fixed cells. 



6 



B 



< 
to 



o 



to 



























_j — 
g 


1 1 1 

20 «0 90 


so too 



Time (sec) 



o 

U 



o 

u 



Fluor reverted/ 
Fluor lost 



1.7S 

1.70 

™ 1 65 

1.S0 

" ' 1.SS 
V 

U 1.58 

V 1-4S 

m 1-40 

V 

US 

o 

zj 1.30 

1.25 

"g 

ra lis 
tjj 1.10 

£i 1 OS 

— Ill 

e.95 




- too 

- 100 

- 100 



20 30 

Ti Time (sec) 



I ■ * :| 



I ■ 



U 



Tau=0.52,sd=0.10 



Tau=34,085, sd= 15,158 



Laser intensity {%) 



Laser intensity (%) 



Laser intensity (%) 



Figure SM3. Reversibility of the photobleaching of YFP when excited with a 515 nm laser. FRAP assays 
where CHX pre-treated cells were completely photobleached. A, example of a reversion of 
photobleaching curve. The blue line represents the "lost" fluorescence and the green, the "reverted". B, 
Curves of reversion of photobleaching for different laser intensities. C, recovered fraction as a function of 
laser intensity. The data was fitted by an exponential curve with parameters: lambda=54,036, A=0,672 y 
B= 0,0548. D, Characteristic times (tau) of the reversion as a function of laser intensity. 



Fluorescence quantification 

To quantify total fluorescence in each compartment, we manually drew regions of interest (ROIs) in the 
nucleus and the cytosol and measured average fluorescence using the ImageJ plugin "Time Series 
Analyzer V2". A typical ROI was a circle of 200 nm in radius (Figure SM4). We studied the effect that 
drawing slightly different ROIs (different shapes, sizes and locations) would have on the estimated 
parameters, but we did not find any significant differences (not shown). 
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Figure SM4. Example of a mother-daughter yeast pair of cells with the ROIs measured. Fluorescence 
corresponds to YFP-Ace2. 



4. Correction for light out of focus 

Because of the optical slicing of confocal microscopes, not all of the light emitted from fluorophores in 
the cell reach the detector when a single slice of the cell is imaged. This effect increases with the size of 
the cell, as a bigger region of the cell is out of focus for bigger cells. The idea is have a measure 
proportional to the total fluorescence of the cell, independently of the cell size. For this purpose we 
simulated images of spherical cells of different size but equal total fluorescence (homogeneously 
distributed). 

To be able to simulate the images, we needed to determine the microscopes PSF. To do this, we acquired 
z-scans of 44nm yellow-green fluospheres (molecular probes). The pixel size is 0.02um in XY and 0.2um 
in Z, laser 515 nm, 12.5 jus per pixel. Images of 512 x 512 px, and 15 Z slices. Because the background 
depends strongly with Z, we acquired an equivalent stack with no sphere. 

We then subtracted the mean background from the image. If lm[i,j,k] is the intensity of the slice k at 
pixel {i,j}, and Bg [i,j,k] is the background image then 

Iml[i,j,k]=Im [i, j, k] - mean (Bg[., ., k]] 

We then corrected for the imaging photobleaching (PB). 

Im2[i,j,k] = e _ (k _ 1)kpb 

where Kpb is the PB rate per image. This rate was calculated by doing an exponential fit to the decay in 
fluorescence of a 

fluosphere imaged always at the same focus. 



To each slice of the Z stack we independently fitted a 2D Gaussian function of the form 

-fcr-xo) 2 -iy-yo) 2 
y] = A * e 2ffx e 2 "y 

This resulted in a table of 15 registers (one per slice) and columns xO, yO, sx, sy, A. Because sx and sy 
behave in a similar way, we averaged both. 
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We tried to fit the PSF to a 3D Gaussian, but failed. The main reason is that s[z] is not constant, and the 
dependence of the amplitude of the 2D Gaussian with z (A[z]) doesn't seem to be Gaussian. We then 
used arbitrary functions to fit s[z] and A[z]. 



We then used the fitted PSF to simulate the light we would acquire by imaging cells of different sizes 
(Figure SM5). When the fluorescence gets recruited in the nucleus, it goes from a occupying a sphere of 
r=2.5mm (the cell) to a sphere of r=lmm (the nucleus). Thus, it gets concentrated in the focal region 
(assuming the nucleus is in focus) and the total fluorescence increases. If we want a measure of the total 
fluorescence, independent of nuclear fraction, we need to multiply to cytoplasm fluorescence by a 
"correction factor". 

The correction factor will be defined as 

^ i-. Fl [Rcell = 1 um] „ „», 
CorrectwnF actor = — : = 1.806 

Fl [Rcell = 2.5 lira] 

In summary, to calculate total nuclear fluorescence, we assumed that we are imaging the whole nucleus, 
which is reasonable approximation considering that using our measured PSF, we established that each 
confocal slice is approximately 1 um, similar to the nuclear diameter. For the cytosol calculation, we 
simulated the microscope optics based on the measured PSF of our microscope (Fig SM5). Using this 
simulation, we estimated that the total fluorescence of the cytosol (for a cell of a 2.5 um radius) is 1.8 
times the signal detected with our microscope settings. 
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Figure SM5. Simulations of imaged cells based on the PSF determined for our microscope. A, Intensity at 
the center of the cell for different cell sizes. As the cell gets bigger the fluorescence gets diluted, and 
more fluorescence goes to out-of-focus regions. B, Integrated fluorescence over the cell area, vs cell size. 
The integrated fluorescence also decreases with the cell size because of the out-of-focus effect. The red 
line shows an exponential fit to the simulation. 



5. Fitting of the nuclear transport model to the FRAP data 

The 2-compartment model 

We used the same model structure, a two-compartment model for all the fittings in the individual cells. 
This model is given by 

d 

— n = —k-i * n + k 7 * c 
dt 

d 

— c = /ci * n — k 7 * c 
dt 
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where n and c correspond to the amount of YFP in the nucleus and in the cytosol, respectively, and ki 
and k 2 are the import and export rates. 

Parameter estimation and optimization 

The experimental data consisted of four sequential FRAPs in the nucleus resulting in four signal 
recoveries. Because of this, the model had to be simulated four times using the same kinetic parameter 
set for the four recovery curves. We formulated and simulated the model in MATLAB using the function 
SBPDsimulate, which can be found in the systems biology toolbox ( http://www.sbtoolbox2.org ). 

In these four different simulations, the initial conditions (IC), i.e. the values that the two states n and c 
starts on, also differ from each other. The IC's in all four simulations were optimized using the first 
observation in the experimental data directly after the perturbation as follows: 

InitialC ' ondition optimized = k IC * InitialC ondition Exver i menta i 

where k IC is the parameter to be optimized. Using the same parameter for the IC in both nucleus and in 
the cytosol was sufficient to describe the experimental data. In total, the fitting of the model output to 
the experimental data from one cell required 6 parameters, 2 kinetic parameters and 4 IC-parameters. 
These 6 parameters where all estimated at the same time. 

For the estimation step, we used the function Isqnonlin, which can be found in the optimization toolbox 
(http://www.mathworks.se/products/optimization/) in MATLAB. Lsqnonlin solves nonlinear least- 
squares (nonlinear data-fitting) problems, and require at least two input arguments; a function FUN, 
whose sum of square is to be minimized and a vector xO, which contains the start guesses of all the 
parameters. Each cell was considered as an isolated estimation problem. 

For this work, we wrote a function FL//Vthat simulated the model for every new parameter set in every 
iteration. Using Isqnonlin, the sum of squares needs not be formed and returned from FUN, since this is 
done in Isqnonlin. Instead, FUN should return a vector containing only the residuals, i.e. the difference 
between the model output and the experimental data, weighted by the experimentally determined 
standard deviation. 

Lsqnonlin is a local optimization algorithm, and because of this, we provided several start guesses of the 
parameters for each cell in order to make sure that the algorithm found the global optimal parameters. 
The alternative would have been to use a global optimization algorithm such as simulated annealing. 
However, using simulated annealing turned out to be more time consuming as compared to having 
several start guesses using Isqnonlin. The start guesses of the IC-parameters were always set to 1, while 
the start guesses for the kinetic parameters where set using a two-dimensional grid ranging from 0 to 10, 
using a step length of 0.1 for both kinetic parameters. After iterating through all start guesses, we chose 
the optimal parameter set based upon the smallest least-square value. Finally, a qualitative analysis was 
made by simulate the model with the optimal parameters and simply study the fitting of the model 
output with the experimental data, which showed a good agreement with all the cells. 
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